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Abstract 

We consider the abelian sandpile model and the uniform spanning unicycle on 
random planar maps. We show that the sandpile density converges to 5/2 as the maps 
get large. For the spanning unicycle, we show that the length and area of the cycle 
converges to the hitting time and location of a simple random walk in the first quadrant. 
The calculations use the “hamburger-cheeseburger” construction of Fortuin-Kasteleyn 
random cluster configurations on random planar maps. 

Keywords: hamburger-cheeseburger bijection, random planar map, abelian sandpile 
model, cycle-rooted spanning tree. 


1 Introduction 


Random planar maps together with discrete statistical mechanics models (e.g., spanning 
tree, Ising model, FK model) on them is an active research area (see e.g., |Ang03 ILGM12I 
iDSm IBenl 31 IShAHl IG.TSZ.T1 21 0041 IDMS141IGMS1 51 IBLR1 51 IChel 51 IGS1 Sal IC1S1 5b] h In 


the annealed distribution of a discrete model on a random planar map of a given class, the 
joint distribution of the pair (M, £), where M is the planar map and £ is a configuration of 
the discrete model on M, is just proportional to the weight of £. Equivalently, the random 
map in the class is sampled according to the partition function of the discrete model (i.e., 
the weighted sum of all the configurations on the given map), and then a configuration is 
sampled according to the weighting rule of the discrete model on the map. 

We consider two discrete models on random planar maps: the uniform recurrent sandpile 
and the uniform spanning unicycle (also known as a cycle-rooted spanning tree). For the 
sandpile model we show that the sandpile density converges to 5/2 and concentrates around 
this value. For the unicycle model, we compute the weak limit of the joint distribution of the 
length of the cycle and the area inside the cycle. The sandpile model calculations depend 
on the results for the spanning unicycle. 


1.1 Planar maps, spanning trees, and unicycles 

Planar graphs are graphs that can be embedded into the sphere. We allow the existence 
of self loops and multiple edges. A planar map is a connected planar graph embedded on 
the sphere considered up to isotopic deformation of the edges, i.e., a planar map contains 
only information about the combinatorial structure of the embedding. A rooted planar map 
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(M, e ) is a planar map M with a distinguished directed edge e. We let M. n denote the set 
of rooted planar maps with n edges. 

Given a finite connected graph Q = (V,E), a spanning forest of Q is a subgraph whose 
vertex set is V and which contains no cycles. A spanning forest that contains k connected 
components is called a k-component spanning forest. A 1-component spanning forest is called 
a spanning tree. A k-excess subgraph of Q is the union of a spanning tree and k extra edges 
in E. A 1-excess subgraph is called a spanning unicycle. The planar dual of a fc-excess 
subgraph is a (k + l)-component spanning forest. Let T(G) denote 

T(G) ■— {spanning trees of Gj ■ 


and £4 (G) denote 

Uk(G) ■— {^-excess subgraphs of Q } . 

Given a rooted planar map (M, e), since e is a directed edge, there is a unique face / to 
the right of e, which we call the outer face. The outer face allows us to distinguish between 
the two sides of a cycle: the outside is the side containing the outer face, and the inside is 
the side which does not. The length of a cycle is the number of edges on the cycle. We define 
the area of a cycle to be twice the number of edges inside the cycle plus its length. Planar 
maps are in natural bijective correspondence with quadrangulations, and this combinatorial 
definition of area corresponds to assigning each quadrangle area 2. 

1.2 The abelian sandpile model 

The abelian sandpile model is a model for self-organized criticality [ BTW88] which is defined 
as follows. (See also [HLMPPW08] for further background.) Suppose G = (V, E) is a finite 
connected undirected graph with loops and multiple edges allowed. Let c(v, w) be the number 
of edges between vertices v and w, where self-loops count twice. For v E V, the degree of v 
is denoted by deg(u) = Y2 w ev c(v, w). A sandpile configuration on the graph G , with respect 
to a distinguished vertex s called the sink, assigns a non-negative integer number of grains 
of sand to each vertex other than the sink s. If a vertex v ^ s has more sand than its degree, 
then v is unstable , and may topple , sending one grain of sand to each neighbor. The sink s 
never topples. Since every vertex is connected to the sink, we may repeatedly topple unstable 
vertices until every vertex is stable. The resulting sandpile is called the stabilization of the 
original sandpile, and is independent of the order in which vertices are toppled (which is the 
abelian property). 

Some sandpile configurations are recurrent, meaning that from any sandpile configura¬ 
tion, it is possible to add some amount of sand to the vertices and stabilize to obtain the 
given configuration. These sandpile configurations are the recurrent states of the Markov 
chain which at each step adds a grain of sand to a random vertex and then stabilizes the 
configuration. The stationary distribution of this Markov chain is the uniform distribution 
on recurrent sandpile configurations. 

We let 7 Z(G, s) denote the set of the recurrent sandpile configurations on a graph G 
with sink s. Majumdar and Dhar gave a bijection between TZ(G,s) and T{G) [ MD91 j. In 
particular, \1Z(G,s)\ = \T(G)\- 

Given a recurrent sandpile configuration o E 7 Z(G, s), where cr(v) for v ^ s is the number 
of grains at v, it is convenient to define cr(s) = deg(s). Then the total amount of sand (i.e., 


2 








including the sand at the sink) is 


M = E a ( u )- 

v€G 

With this convention, the distribution of |cr| for a random recurrent sandpile a G lZ(G,s) 
does not depend on the choice of the sink s |Mer97j . This distribution can be understood in 
terms of the k-ex cess subgraphs of G |Mer97j , which we will explain in Section [2J 
We define the sandpile edge density for Q to be 


Pe(G) 


1 

E\-\R(G,s) | 


E 

cr£'R.(G,s) 


a 


and the sandpile vertex density for Q to be 


Pv(G) 


1 

WF\mW\ 


E 

&£7l(G,s) 


a 


( 1 ) 


( 2 ) 


These sandpile densities p e {G) and p v {G) are independent of the choice of the sink. 


1.3 Sandpiles and k-ex cess subgraphs on random planar maps 

We let denote the annealed distribution on ^-excess subgraphs of a rooted planar map 
containing n edges, i.e., a random planar map (M, e) with a k-excess subgraph 74 on it are 
chosen with probability 


KlM,e,U h] = E i - wm . (3) 

{M',e')eM n 

We let E(j and Var^ denote the expectation and variance with respect to P(j. 

For the directed edge e, we let e denote its source vertex, and e denote its destination 
vertex. Because of the bijection between recurrent sandpile configurations 7 Z(M,e) and 
spanning trees T(M) = Uq (M), we can interpret as being the annealed distribution of 
uniform recurrent sandpiles on A4 n with sink at the source of the root edge. 

Theorem 1.1. For the uniform recurrent sandpile on a randoin planar map with n edges, 
the sandpile density satisfies 

hm E“[p,,] = 2 lim E ° n [p e ] = lim Var ° n [p v \ = lim Var°[p e ] = 0. (4) 

This sandpile density computation can be compared to the uniform recurrent sandpile 
density on Z 2 . In 1994 Grassberger conjectured that the (per vertex) sandpile density on 
Z 2 is 17/8, based on the numerical integration of singular 4-dimensional integral expressions 
given by Priezzhev |Pri93j for the sandpile height distribution at a vertex. These integral 
expressions were greatly simplified by Jeng, Piroux, and Ruelle [JPR06J, who verified by 
numerical integration that the sandpile density for Z 2 is 17/8 ± 10~ 12 . An alternative for¬ 
mulation of the sandpile density, relating it to spanning unicyclic graphs and loop-erased 
random walk, was given by Poghosyan and Priezzhev [PPllj and Levine and Peres jLP 1 4] , 
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which enabled its rigorous exact evaluation in [PPRlll IKW14aj . The density of 17/8 corre¬ 
sponds to an edge density of 1 + ^, versus 1 + f for the sandpilc on a random planar map. 
More recently, a simpler proof of the sandpile density on Z 2 was given by Kassel and Wilson 
|KW14bj . who then computed the sandpilc density for numerous other lattices as well. 


Our results for unicycles, and more generally k-ex cess subgraphs, can be expressed in 
terms of simple random walk in the first quadrant of Z 2 . Suppose (X t ,Y t ) is a simple 
random walk on Z 2 started from (1,1). Denote 


^quad — inf{7 : X t Y t — 0}, s quad — ^D quad + Pt quad • (5) 


In words, f qua d is the exit time from the first quadrant, and s qua d is the distance from the 
origin when the walk exits the quadrant. Random walk in the first quadrant has been studied 
quite extensively (see e.g., ( FIM99 . Rasl2 . DW15 . IRLIOj ): for an introductory account see 
(LL101 Section 8.1.3]. Here we record that 


P[*quad > j] ~ 


4 

Tlj 


and P[s quad > £] ~ 


( 6 ) 


asymptotically as j —* oo and £ —* oo (see [ FR12 


Sec. 3] and [DW151 Example 3, Sec. 1.3]). 


Theorem 1.2. Consider the uniform k-excess subgraph of a random planar map under the 
distribution ofF (j. As n —» oo with k fixed, with probability 1 — o n (l) the k-excess subgraph 
has k unnested loops. Let L\, Ai ,..., L*,, denote the length and area of the k loops. Then 
the joint distribution of {(Li, Aj)}i<j<fc converges weakly to k i.i.d. samples of the random 
variables (s qua d, t qua d) described above. 


The area distribution can be compared to the area of the cycle of a uniform spanning 
unicycle in Z 2 . The moments of the area were computed for n x n boxes in Z 2 [ KKW12] , 
and these moments suggest that IP [A > j} x 1/j. 


To prove these theorems, we start in Section [2] by explaining how the k th moments of 
the amount of sand in a uniform recurrent sandpilc are related to k-excess subgraphs. Then 
in Section [3] we review the “hamburger-cheeseburger” bijection, which constructs rooted 
planar maps with n edges together with a Fortuin-Kasteleyn configuration. (These FK 
configurations are more general than fc-excess subgraphs.) This bijection is due to Mullin 
[Mul67j in the special case where the FK model is a spanning tree, and to Bernardi and 
Sheffield }Ber08l IShellj in general. We use the formulation in (Shell] since it is more 
convenient for our purposes. In Section [4] we do some asymptotic analysis of the hamburger- 
cheeseburger bijection to prove Theorem 11.21 We conclude the proof of Theorem 11.11 in 
Section 0 


Acknowledgments: We thank Adrien Kassel for helpful comments. This work was begun 
while the first author was an intern in the theory group at Microsoft Research Redmond, 
and was completed at the Newton Institute. 
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2 Sandpiles, spanning trees, and the Tutte polynomial 


One approach to computing the sandpilc density is via the Tutte polynomial. The Tutte 
polynomial of an undirected graph Q = (V, E) is a polynomial in two variables defined by 

Tg(x, V) 1 (y - I)“l p >+|E'1-I^| (7) 

E'CE 

where k(E') is the number of connected components of the spanning subgraph of Q with 
edge set E'. 

For the abelian sandpile model on a finite connected graph Q = (V, E) with sink s, Biggs 
defined the level of a sandpilc configuration to be 


level(er) = |cr| — \E 


and showed that 0 < level(cr) < \E\ — \V\ + 1 and that these bounds are tight. Consequently, 
the sandpile edge density satisfies 

1 < Pe < 2 • 

Biggs conjectured and Merino proved |Mer97] that the generating function of recurrent sand- 
piles by level is the Tutte polynomial evaluated at (1,?/): 


E, =Tg(l,y). (8) 

crGTZ(G,s) 


Notice that the generating function is independent of the choice of the sink s. 
Since we are assuming that Q is connected, from (J8J) and (J7J) we see 

/level (a) 


£ l . 

sandpiles a of Q x 

so for a random recurrent sandpile a G 7 Z(G, s ), 

/level(a) 


E e 




l 


= W9) I, 


mo)\ 

\T(G)\ 


(9) 


( 10 ) 


When the graph Q is itself random, in particular a random planar map on n edges, we have 


/level(cr)\ 

= E° 

\\Ui(M n )[\ 

LI * )\ 


[\T(M n )\\ 


( 11 ) 


This equation relates the binomial moments of the sandpilc level to quantities that can be 
evaluated by the hamburger-cheeseburger bijection which we will describe in Section EJ 
For any random variable Z, the moments E \Z l \ of Z can be expressed as linear combina¬ 
tions of the binomial moments E[(^)], and the cumulants E [(Z — E[Z]) £ ] can be expressed 
as linear combinations of products of the form E[Z m ]E[Z] £-m . The variance in particular is 
Var[Z] = 2E[(f)] + E[(^)] — E[(^)] 2 . For a random recurrent sandpile a on a random map, 


E° [level (cr(M n ))] 


Var° [level (a(M n ))] 


, 0 \ mM n )\ - 

,n [\UM n )\_ 

\Ui (M n )\ - 
\T(M n )\ _ 


2 E° 


\U 2 {M n )\ 

I \T(M n )\ I 


+ E 


/ 0 r \Ui (M n ) | ~ 
V n l\T(M n )\_ 


2 


( 12 ) 

(13) 


Next we evaluate the terms on the right-hand side using the hamburger-cheesburger bijection. 
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3 The hamburger-cheeseburger bijection 


Sheffield (Shell! Section 4.1] constructed a bijection called the hamburger-cheeseburger bijec¬ 
tion between “perfect words” over a five-letter alphabet { (H), [R], ©, ©, [F]} and FK configu¬ 
rations on rooted planar maps. The letters in this alphabet can be interpreted as events at 
a burger restaurant. (H) and © indicate that a new hamburger or cheeseburger is produced. 
New burgers are placed on top of a burger stack, which is initially empty. [R] and [C] indicate 
that a customer has ordered a hamburger or cheeseburger, in which case the topmost burger 
of the appropriate type is removed from the stack and given to the customer. |F] indicates a 
fresh order, where the customer orders whichever burger is on top of the stack, regardless of 
type. A perfect word is a sequence of these letters for which every burger order is fulfilled by 
a burger already on the stack, and for which the burger stack ends empty. A perfect word 
of order n is one which contains n burgers and n orders, i.e., has length 2 n. 

The hamburger-cheeseburger bijection maps a perfect word of order n to a rooted planar 
map (M n ,e) with n edges together with a subset E' of the edges in the map. (The edge 
subset E' is the FK configuration on the map.) In addition to Sheffield’s description of the 
bijection (Shell! Section 4.1], a nice exposition is given by Chen (Chel5j . We note here a 
few basic properties of the bijection: 

1. The number of edges in E' is the number of [R]’s plus the number of |F]'s matching 


1 s. 


2. The number of connected components in the subgraph spanned by E' is 1 plus the 
number of [F]’s matching (H)’s. 

3. Let (. E')* denote the dual edges of E\E' on the dual map M* of M n . The number of 
connected components in the dual subgraph spanned by ( E')* is 1 plus the number of 
[FI’s matching ©’s. 

4. Suppose that a |F] matches a © in a perfect word, and that in between there are £ 
Til’s which are fulfilled by (H)’s before the fresh ©, and 2 m other letters. Then the [F] 
corresponds the “last edge” of a loop of E' which has length t + 1, the portion of the 
map M n inside the loop has area (1 + 1 + 2m), and the portion of M n and E' inside of 
the loop are determined by the subword between the |F] and its matching ©. 

Proposition 3.1. Let 

(-)(] = { perfect words of order n with exactly k pPps, which are all fulfilled by ©’s} . (14) 

Under the hamburger-cheeseburger bijection, elements in ©Jj correspond to triples (M n , e, Uk), 
where (M n ,e) € A4 n and Uk is a k-excess connected subgraph of M n . Furthermore, 


El 


\U k (M n )\ 


1 © 

|0 


C \k I 
n I 

01 
n I 


(15) 


V\T(M n ) | 

Proof. Suppose E' has 1 component and its dual (. E')* has k + 1 components. Then (. E ')* 
contains no cycles, so it is a (k + l)-component spanning forest, and consequently E' is a 
k-excess subgraph of M n . 

Recall that under the distribution P°, each rooted map (M n , e) occurs with probability 
|T(Af n )|/|T(M')| = |T(M n )|/|©°|. In the expectation in (JE2D, the |T(M n )| 
terms cancel, giving the right-hand side of (1T5|1 . □ 
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4 Asymptotic enumeration of perfect words 

In this section we compute the asymptotic number of perfect words in |0*|, which together 
with equations (fT2]h (TT3|) . and (fT5li gives the sandpilc edge density. In the course of charac¬ 
terizing 0^, we also characterize the cycles in k-ex cess graphs. 


Theorem 4.1. For any fixed nonnegative integer k, 


lim 

n—>oo 



l 

k\ 4 fc 


(16) 


To prove Theorem 14.11 we study the canonical injection from to 0° x Given 

a perfect word W G 0)j, we can replace each of the k [fj’s in W with ©s to obtain a 
perfect word in 0°, and record the positions of the [F]’s as a set of k distinct elements of 
[2 n] = {1,... , 2n}; this map is invertible so it is an injection. A word W in 0° together 
with a set of distinct positions {A,... , A} is i 11 the inrage of this injection precisely when 
W tl ,..., Wi k are all ©s, and just prior to each of these orders, the top burger in the stack 
is a ©. 

Theorem 14.11 can be interpreted as a statement about the probability that a random 
element of 0 ° x (^) is in the image of the injection. For a random word W G 0 ° and a 
random position i, P[Wj = C ] — and at a random time, provided the stack is nonempty, 
the top burger on the stack is a © with probability |. For large n, it is plausible that these 
events at the same random time are approximately uncorrelated. As long as k is not too 
large ( k <C \/n), we expect random distinct positions A, ■ ■ ■ ,4 f° be far apart, and that 
consequently these events at the times A, • • ■, A are nearly independent. Provided that this 
intution is correct, then |0jj| fa |©n| (©)/8 fc , which when k <C \fn would give the theorem. 
In this section, we justify a more precise version of this intuition to prove the theorem. 

To make a more precise statement of this approximate independence, we consider sub¬ 
words of W G 0°. Let W[a,b] denote the subword from positions a through b inclusive. Let 
Wj be the subword 

Wj := lF[max(A — s + 1,1), A] (17) 

for j = 1,..., k, i.e., Wj is the subword of length s which ends at position ij (unless the 
position is too close to the front, in which case the length will be less than s). Since the 
perfect word W corresponds to random walks in the quadrant that start and end at the 
origin, we expect the subwords wi,...,Wk to be close in distribution to i.i.d. uniformly 
random words of length s, so long as both k and s are small enough for the subwords to be 
disjoint and not to contain enough letters to detect that W is not quite an unbiased random 
walk. 


Theorem 4.2. Assume k 3 s 2 -C n, and consider the collection of subwords (w\,... ,Wk) 
defined in C3 from a random perfect word in 0° and independent random indices A < 

• • • < A in (bG). The subwords are nonoverlapping with probability 1 — o(l) and have total 
variation distance o(l) from a list of k i.i.d. uniformly random words of length s which are 
independent of the indices. (The o(l) terms go to zero as k 3 s 2 /n —> C© 


Proof. Let A < • • • < A be a random /c-tuple uniformly drawn from (^), and let W\,... ,u>k 
be k i.i.d. uniformly random words in {(H), ©, [R], ©} s , which are also independent of the 
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Our goal is to sample a uniformly random perfect word in W G 0° so that 
the subwords defined by (|I71) using the indices ii,...,i k coincide with uq,..., w k - Our 
strategy is to first sample an independent uniformly random perfect word X G 0°, and then 
modify it, using A,..., i k and u>i ,..., w k and some auxiliary randomness, to obtain W. This 
modification proecedure defines a Markov chain, and we restrict ourselves to modification 
procedures for which the Markov chain satisfies detailed balance conditional on A,..., i k , i.e., 
P[A — >■ W | ik] = P[W —>■ X | ii ,..., i k ], so as to ensure that W is uniformly random, 

and in fact uniformly random even conditional on ii,... ,i k . After verifying detailed balance, 
we argue that with high probability W has the desired subwords: w\ = w\,... ,w k = w k - 
We use the following modification rule: If i\ < s or i J+1 — ij < s for some j, then no 
change is made. Otherwise, all the subwords have length s and are disjoint. We then start by 
overwriting the relevant positions in X with the values given by W \,..., w k to obtain a new 
word Y. The walk in Z 2 defined Y is unlikely to return to its start, so some more changes 
are required to rebalance it. 

To describe this rebalancing of the walk, we follow Sheffield [Shell] in using a pair of 
coordinates that are rotated 45° from the edges of Z 2 : we let u denote the number of 
burgers minus the number of orders, and v denote the “discrepancy” between hamburgers 
and cheeseburgers: 

#{©,©} - #{0,©} = u and #{©,B} - #{0, Q} = v. 


In these coordinates, the letters correspond to the following steps: 


(ED — ( + 1) +1) © — ( + 1; — 1) 


0 = (—i, +i) 


H 


(- 1 ,- 1 )- 


If u > 0, for instance, then focusing on the first coordinate, we change u/2 of the +l’s to 
— l’s, while ignoring the second coordinate. The second coordinate can then be rebalanced 
ignoring the first coordinate. When doing this rebalancing, we only change letters whose 
position is in the range from to | n, and which do not lie within the subwords. For reasons 
that will become apparent, out of all possible such ways to rebalance the walk Y, we pick 
one uniformly at random, and let Z denote the resulting walk. (If there are no such ways to 
rebalance Y, we let W — Z — X.) 

We will argue later that Z is likely to remain within the quadrant, but if not, then we 
let W = X. Otherwise, Z is a perfect word, and we would like to take W — Z, but to 
ensure detailed balance conditional on i±,..., ik, we consider Z to be a proposal, and use the 
Metropolis rule to reject this proposal with some probability and instead take W = X. The 
probabilities that X proposes Z and that Z proposes A" (conditional on A,..., ik) are almost 
the same, but there is a small difference arising from the (likely) possibility that there are 
different numbers of ways to do the rebalancing in the two cases. If there are Tq ways to do 
the rebalancing when going from X to Z, and r 2 ways when going from Z to A", then the 
probability of accepting the proposed move is min(l, ri/r 2 ). 

We have specified the process which produces the uniformly random perfect word W, 
and at this point we argue that with high probability all the steps succeed so that in the end 
w 1 =w l ,...,w k = w k . 

Since k 2 s <C n, with probability 1 — o(l) the subwords are all disjoint and have length s. 

After the positions are overwritten to obtain Y, certainly |u|, |u| < 2 ks. Since ks -C n, 
the number m of letters with position between |n and |n outside the subwords is certainly 





m = (l + o(l))n. For each of the two coordinates, standard large deviation estimates together 
with the cycle-lemma construction of random Catalan paths [DM47 ] imply that there are 
between | m — ay/m and + ay/m letters with +1 in the first coordinate (and similarly 
for the second coordinate), with probability tending to 1 as a —)■ oo. Assuming this event 
occurs, there is a way to rebalance the walk. Provided that the rebalanced walk Z remains 
in the quadrant, since |w|, |v| < 2 ks yfm, the ratio ri/r 2 tends to 1, so the proposed move 
would be almost always accepted. 

We are left to argue that the rebalanced walk Z almost always remains in the quadrant. 
Suppose that between times en and (2 — e)n the walk X remains at distance at least h from 
the boundary of the quadrant. If £ l/k and s < en, then with probability 1 — o(l) all of 
the subword regions are contained with the interval from en to (2 — e)n. There are at most 
3 ks letters that get changed. If h > 6 ks, then we would be guaranteed that the modified 
walk Z would remain in the quadrant. 

For the initial perfect word X, if we eliminate the © and © letters, will be a random 
Catalan path of a random length 2£ which is concentrated around n ± 0(y/n). Consider a 
uniformly random Catalan path of length 2 L As £ —> oo, near its endpoints it behaves as 
a Bessel(3) process. Since a Bessel(3) process always remains positive, between positions 
e£ and (2 — e)£ the Catalan path is likely to be at height at least 5\fel, with probability 
tending to 1 as 5 —> 0. Setting e = o{l/k) and 5 = o(l), we find that X remains distance 
at least 6 ks from the boundary of the quadrant with high probability provided ks y/n/k, 
i.e., k 3 s 2 -C n. □ 


Using this approximate i.i.d. property of the perfect words in 0°, we can characterize 0(j 
to prove both Theorem 14.11 and Theorem 11.21 


Proof of Theorem\f.l\ and Theorem \1.2\ So long as k 3 s 2 <C n, Theorem 14.21 gives a coupling 


with 


of the subwords W\,...,Wk of W €0° at random locations A <•••<?& G ( ^ ) 
i.i.d. uniformly random words of length s, so that w.h.p. the subwords equal the random 
words. In the i.i.d. words, the probability that the last letter is © is As s gets large, 
the probability that a given subword does not contain a burger that is put on the stack just 
prior to the © is 0(1/y/s), and the probability that this burger is a © is We can extend 
the i.i.d. words by prepending uniformly random letters, and then the probability that for 
each word the last letter is a © that matches a © is exactly 8~ k , and with probability 
1 — 0(k/\/s ) each matching © occurs within the (unextended) subword. With probability 
8~ k — 0{k/ ©s) the word W with the inidices (A, ..., A) is in the image of the injection, and 
provided k/ ©s 8~ k , this is 8 -fc (l — o(l)). This proves Theorem 14.11 

Assuming again that k 3 s 2 n and k/y/s -C 8~ k , the above coupling gives a character¬ 
ization for a random perfect word W k in 0(j. Suppose the fresh orders occur at positions 
i\ < ■■■ < ik, and ji,---,jk are the positions of the corresponding fresh cheeseburgers. 
Then with high probability these indices alternate, j\ < A < j 2 < A < ■ • • < jk < h, so 
that the loops given by the hamburger-cheesburger bijection are unnested, and the subwords 
W k \j\-\-l, A — 1], • • •, W k [jk T1, A — 1] are within o(l) variation distance of k i.i.d. simple ran¬ 
dom walks of the following type: when reading backwards, the hamburger and cheeseburger 
counts are always nonpositive, and the cheeseburger count ends at 0. These walks are of 
course equivalent to walks in the quadrant, and from the discussion in Section [3l the length 
of the walk corresponds to the area of the cycle, and the hamburger deficit corresponds to 
the length of the cycle. This proves Theorem 11.21 □ 
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5 Sandpile density 


The sandpile (edge) density is a straightfoward consequence of the above lemmas: 

Proof of Theorem M. 1\ (edge density). Combining equations (TT2]h (fT3]h (TT5]h and (1T6l) gives 
E° [level(cr(M n ))] = (1 + o(l))n/4 and Var° [level(cr(M n ))] = o[n 2 ). The amount of sand 
is n more than the level, which gives the sandpile density with respect to the number of 
edges. □ 

For the random planar map with n edges, the expected number of vertices is n/2 + 1. 
To obtain the sandpile density with respect to the number of vertices, we need to know that 
the number of vertices is sharply concentrated about its expected value. 

Lemma 5.1. For any fixed k > 0, for a random planar map M n drawn from (A4 n , P°), 


lim E° 

n—>• oo 


fjmoi 

\2| V (A/„)| 


k 


l. 


Proof. Let J denote the number of (H)’s contained in a uniformly random perfect word in 
0°. With Catj denoting the jth Catalan number, we have 


|0°|xP[J = j] = 


2 n 
2j 
n + 1 
3 


CatjCat n _j = ------ 

j\(j+l)\(n-j)\(n-j+ 1)\ 


n + 1 
3 + 1 


(2 n)\ 


(n + l)!(n + 1)! 


As is well known, the binomial coefficients are sharply concentrated, with tails that are 
at least as small as Gaussian tails. In particular, there are positive constants C and c for 
which 

P[|J-a±!| >t]< Ce~ ct2/n . 

Recalling the properties of the hamburger-cheeseburger bijection from Section [3j the 
number of edges in E' is | ,/|, and since E' forms a spanning tree, the map M n has \V\ = | J\ + 1 
vertices, while of course the number of edges is \E\ = n. Regardless of how atypical | J\ may 
be, we have 1 < \V\ < n + 1, so n/(n + 1) < |£ , |/|R| < n. With, say t = n 2 / 3 , we obtain 


E [(\E\/\V\) k ] = E[(|E|/|R|) fc l ||F| _ n/2|> ,] + E[(|£?|/|y|) fc l|| V| -n /2 |< 
= n fc e“ 0(nl/3) + 2 fc (l + 0(kn~ 1/3 )). 


□ 


Proof of Theorem \El\ (vertex density). It suffices to prove that lim, woo E°[(p„ — 5/2) 2 ] = 0. 


Pv 


\= 2 P . 


E I 


Gvi 

P ’ : ~ 2 ) _4p '(w 


+ 2 ' Pe ~ 4 


+ 8 p e \p e -~ 


\E\ 


1 I + 4 [ Pe~ ~ 


Recall that 1 < p e <2. By Lemma Hull lim^,*, E° [(|FJ|/|R| — 2) k ] = 0 for any fixed integer 
k > 0, so the first two terms above converge to 0 in expectation. From the edge density part 
of Theorem 11.11 we see that the last term converges to 0 as well. □ 
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